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8 : Abstract 



The lepton number violating signal of like-sign-dileptons (LSD), pp — ► £ £ + 

o 

m 

C the existence of heavy Major ana neutrino states with a left-handed coupling 



2jets, is investigated within a model of mixing in the neutrino sector assuming 



to the light leptons. The LSD signal receives contributions both from the 

resonant production of a heavy Majorana neutrino (N) and from the exchange 

P-t. 

D . of a virtual iV in the WW fusion mechanisms. These two possibilities are 

discussed in detail and compared. Helicity amplitudes are given pointing out 

differences with calculations previously reported by other authors. The signal 

cross-sections are computed at the energy of the LHC collider (v5 = 14 TeV) 

at CERN and within the existing experimental limits on the mixing couplings, 

including those coming from neutrinoless double beta decay. Detailed angular 

distributions of signal reactions which are complementary to previous studies 

on the argument are presented. 
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I. INTRODUCTION 

The quest for a complete understanding of the true nature of neutrino masses and mixing 
in the neutrino sector of the Standard Model (SM) is one of the most important challenges 
that physicists have to face in the coming years. Perhaps the experimental difficulties 
involved in the study of neutrino properties can be best appreciated noticing that only 
very recently the neutral partner of the tau lepton, the tau neutrino (v T ), nas been directly 
observed [IJ. In addition, after decades of experimental efforts there is finally some evidence 
for neutrino oscillations and thus the first indications of the fact that neutrinos might 
not be massless particles as advocated by the SM. Even so, it is not clear at present whether 
these neutral fermions are of Dirac or Majorana type (for a Dirac neutral v 7^ v while for 
a Majorana v — v). It is well known that if neutrinos are of the Majorana type then the 
Lepton number (L) is not conserved and processes which violate L by two units AL = ±2 
become possible. Perhaps one of the most interesting lepton number violating processes is 
the neutrinoless double beta decay j3(3q v which currently provides the most stringent bound 
on the effective light Majorana mass (m v ) < 0.1 eV. Extensive literature is available on this 
subject both from the experimental and theoretical sides. In view of the forthcoming large 
hadron collider (LHC) to be built at CERN which will realize proton-proton collisions at 
energies of v5 = 14 TeV it will be of increasing importance to study quark reactions which 
are the high-energy analogues of the /3/3o u - This was realized already some time ago and 
has been the object of several studies of the signal (lepton number violating) reaction: 

pp -> &t + 2 jets (AL = ±2). (1) 

Any model with lepton number violation (i.e. due to the presence of massive Majorana 
neutrinos) will provide a non zero amplitude for the process in Eq. [I], which receives con- 
tributions from parton scattering as depicted in Fig. p] and Fig. ^|. In a series of recent 
papers 0-0, the present authors have addressed the question of lepton number violation 
within a composite model scenario where the partner of the light neutrino, the so called 



excited neutrino, N, is assumed to be of a Majorana type. The phenomenological impli- 
cations of this idea have been investigated in detail both with respect to the low-energy 
/3/3ov decay and to the high-energy LSD signal in hadron collisions. In particular, in || we 
developed a helicity amplitude method to efficiently compute all 2 — ► 4 parton subprocesses 
that contribute to the reaction in Eq. ([!]). The cross-sections were estimated to lead, within 
some regions of the parameter space, to an observable signal. One particular advantage 
of the approach developed in || is that it gives the possibility to treat coherently the two 
mechanisms of WW-fusion and resonant production (annihilation channel). Clearly such 
an approach can be extended to explore the phenomenology of other models that entail L 
violation. In the event of observation of such a signal, it will be important to know the 
contribution to it of the various models beyond the SM and ultimately, a way to distinguish 
among them to understand the responsible mechanism . 

In this work, the model of heavy isosinglets Majorana states that mix with ordinary light 
lepton via a left-handed coupling suppressed by a mixing coefficient, is assumed and studied 
in detail. While completing this work a report appeared in the literature which also 
commented upon the same model of heavy Majorana neutrinos by studying the resonant 
production at LHC via the third order process 

pp -> £ + £ + + W~ AL = +2 (2) 

and arriving at somewhat optimistic results as to the possibility of observing the signal. 
The process in Eq. (|j) is described by the Feynman diagram obtained by that of Fig. |] but 
without including the hadronic decay products of the final on-shell W. However, the method 
used here of evaluating analytically all the relevant amplitudes, as developed in ||, is totally 
independent of that of the authors of ]/]] which relies on the CompHEP || package. The 
rest of the paper is organized as follows: section II gives a description of the theoretical 
framework; in section III the amplitudes of the contributing subprocesses are presented; 
in section IV an exact analytic result is provided for the process in Eq. [2] highlighting 
differences with the formula reported in [0] and discussing thoroughly how the conclusions 



reported there are affected; section V presents some discussion of kinematic cuts and gives 
some details of the numerical calculations along with the results; finally section VI provides 
concluding remarks. 

II. THE MODEL 



The mixing of heavy and light neutrinos is a common feature of many theoretical models 
beyond the SM. Heavy Majorana neutral leptons have been well studied in the literature in 
particular from a model which is the SM augmented with additional right-handed neutrino 
fields. In this scenario, that can be seen as a low-energy manifestation of SU(10) GUT 
theories P,|T0||, one has the following mass term for neutrinos: 
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where {vl) c and vr are vectors in generation space with components given by three chiral 
fields: 
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and a mass matrix given by 
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In general an arbitrary number of right-handed neutrinos can be added, but for definite- 
ness and simplicity the most symmetrical case is considered (i.e. that of three right-handed 
fields). The mass matrix is then 6x6, with Md the 3x3 matrix of Dirac type masses for 
neutrinos that is built via the spontaneous symmetry breaking and Yukawa couplings as for 
charged leptons and quarks, while Mr is the Majorana 3x3 mass matrix for singlet right- 
handed fields. If one requires the eigenvalues of Mr to be much bigger than those of Md, 



diagonalizing the mass term one has a three generation "see-saw" ||il| , p!2| model, with three 
light and three heavy Majorana mass eigenstates. Weak interactions produce weak eigen- 
states that are superpositions of mass eigenstates through the mixing matrix. As opposed 
to the naive one family "see-saw" mechanism, where 9 ~ \Jm v / M^, with three generations 
(or more) of right-handed neutrino fields, it is possible to decouple the magnitude of the 
mixing coefficient of heavy-light states from the mass eigenvalues of Mr, and the mixing 
angles can be treated as phenomenological parameters bounded only by existing experimen- 
tal data. Examples of mass matrices that satisfy this request are given in ref. [[13,141,9,15 



In this context, lepton number violation in pp collisions was analyzed in |^,|l3|,[n]] . Some 
authors refer instead to other theoretical scenarios [f7|)|18| derived from super-string inspired 
Eq models. Here global symmetries imposed on the mass matrix decouple the mixing angles 
from any relation with masses eigenvalues, so that they are free parameters that can be large 
(as in multi-generational "see-saw" type models). Some caution however must be paid in 
deriving heavy Majorana masses in these models where every family of the SM is enlarged 
both with a right-handed field and a singlet field Si in the neutrino sector. The mass term 



for neutrinos in the Lagrangian is [19-21 
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where the mass matrix is given by: 
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where again the fields represent a collection of three chiral fields and the block M and D are 
Dirac type sub-matrices. The zeros are consequence of the global lepton number conservation 
imposed on the mass matrix. The spectrum that one obtains through diagonalization is 



three massless neutrinos and three heavy Dirac neutral leptons p2 |. Therefore these type 



of models cannot be the base for the discussion of the lepton violating processes that are 
the object of this work. One could add a Majorana mass matrix, /i for Sl [[£J (lower left 
3x3 block in the mass matrix of Eq. (|7])), and obtain a small Majorana mass for the light 
neutrinos, but the heavy states remain Dirac particles. Thus in this class of models it is not 
very easy to obtain heavy Majorana masses (in the TeV range) with substantial heavy-light 
mixing. 

Independent of theoretical prejudices, in the following the general interaction Lagrangian 
between charged leptons, neutrinos, and gauge vector bosons of the SM is given within this 
class of models where the weak eigenstate neutrinos are combinations through the mixing 
matrix of light and heavy physical neutrinos. So the latter have the usual left-handed 
coupling with leptons with a coupling constant that is the SM weak coupling g multiplied by 
the appropriate element of the mixing matrix. The charged and neutral current interaction 
Lagrangian densities are given by: 
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where Bi n is a 3 x 6 rectangular mixing matrix that can be written Bi n = (£?&,, Bin), with 



Bi v and B^ the two 3x3 block matrices that give the light-light and light-heavy mixing, 
while the 6x6 matrix CV, is given by C = B^B. The properties of B and C are extensively 



discussed in the literature [24 -661 . We now discuss the known experimental constraints on 



mixing angles. As these heavy states have not yet been observed, light-heavy mixing has to 
be inferred by low-energy phenomenology. In [27], a global fit is performed on LEP data 
identifying the following effective mixing angles: 

c 2 Ve = {bIb 1v ) u = J2 \{B iv ) u \ 2 =cos 2 e ve , 
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that are bounded from above as: 

si < 0.005 si < 0.002 si < 0.01. (10) 

The authors of ref. have performed an updated analysis finding roughly the same value 
for the bound on the electron parameter s 2 Ue < 0.0052 but a somewhat more stringent 
bound on the muon mixing parameter s 2 u < 0.0001, thereby making the electron LSD 
channel the most interesting. However for first generation leptons it is necessary to consider 
bounds coming from non observation of neutrinoless double beta decay (see below). In the 
next sections calculations are reported for the first family (£ = e) LSD signal under the 
assumption that: i) only one heavy Majorana eigenstate is exchanged; ii) the maximum 
allowed value for the mixing can be used for the relative element of the mixing matrix, i.e. 
|-Be7v| 2 = 0.0052. It should be recalled that when extensions of the SM are embedded in more 
general theories, one needs to consider new gauge bosons that can interact with the heavy 
neutrino states. However, as these bosons are predicted to be very heavy, their interactions 
are suppressed relative to those of the SM so that one can safely neglect them |7|,[| in the 
following discussion. 

A. Unitarity and (3(3q u constraints 

Before entering into the detailed calculations, a comment is in order, regarding violation 
of unitarity, and the bound coming from the non observation of neutrinoless double beta 
decay (/3/3 „) which is quite important. It is well known |T6| , |28|JT8| that the basic process 
e^e - — > W~W~ or, equivalently, W~W~ —> e~e~ intervening in the reactions discussed in 
the present work (see Fig. [l]), violates unitarity unless the following condition is satisfied: 

Y,m v {B iv f + Y. Mn ( B ^ 2 = ° ' ( n ) 

v N 



In the model considered here this is automatically satisfied as it is easily shown that 
~Yli V m v (Bi u ) 2 + ^2 N M N (B iN ) 2 = M u where M K is the left-Major ana mass matrix (3x3 
block matrix) of vg, and this is zero as can be inferred from the mass matrix given in Eq. (0). 
The Majorana mass term of U£ could be different from zero within a Higgs triplet model 
which however has been excluded by the precision LEP data, unless one considers somewhat 
more contrived models, where, for example, the Higgs triplet is mixed with a singlet allowing 
to evade the LEP bounds p9| . In the following such eventualities are not considered, the 
absence of a Higgs triplet is assumed, ensuring thus the unitarity condition to be satisfied 
according to Eq. (H). 

The non observation of the neutrinoless double beta decay (a second order weak nuclear 
decay) gives the following bound on the effective inverse mass |28|| : 
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< 7x 1(T 5 TeV" 1 , (12) 



Within the rather reasonable hypothesis that only the light mass eigenstate contributes 
significantly to both the low-energy (P/3q u ) process and high-energy (LSD) signal and using 



|S e jv| — 0.0052 one easily finds a rather impressive lower bound on the mass [28 



M 7V >75TeV. (13) 

This of course would preclude the possibility of producing such particles at any forthcom- 
ing hadron or lepton collider. The (3f3 0v bound could obviously be evaded as the mixing 
coefficients appear as squared and not through their moduli squared in Eq (|T2"D. Indeed 



in ref. [28| several scenarios of mixing are discussed where cancellations between differ- 



ent terms arise thus evading the PPq v bound. Explicit examples are provided [T2I| but of 
course the resulting models appear somewhat fine tuned and/or too much contrived and 
un-natural. It would then appear that a phenomenological analysis of the processes dis- 
cussed in the introduction would be pointless for heavy neutrino masses within the mass 
range 100 GeV < Mn < 1 TeV. However, one can build models which do not suffer such 
drawbacks of fine tuning and unnaturalness by simply assuming that the heavy Majorana 
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neutrinos do not have the same CP parities (rjcp)- This was done by the authors of ref. [18 
who considered a model with two degenerate heavy Majorana neutrinos of Mjy ~ 1 TeV 
and opposite CP parities that evades the (3(3q v bound |18[. It is within such a scenario 



that the following phenomenological study of the processes in Eq. (|2]) and Eq. (|1|) has to be 
understood. 

In any case it must be kept in mind that the predictions reported here are based on the 
maximum experimentally allowed mixing and therefore new global fits on the last LEP data 
may lower its value, and thus affect the present conclusions. 

III. AMPLITUDES OF L- VIOLATING PARTON SUBPROCESS 

In the following, the helicity amplitudes for parton sub-processes that contribute to 
production of LSD via the exchange (or production) of a heavy Majorana neutrino are 
presented. The effective interaction used is that of Eq. (|8|). Considering for the moment only 
the first family, three different types of processes should be distinguished: {%) uu — > dd+£ + £ + , 
(ii) ud^ du + £ + £ + , (Hi) dd -> uu + £ + £ + . 

Let the tensor T^ u describe the virtual sub-process W*W* — ► £ + £ + (Fig. [l]), while the 
tensor T^ v describes the virtual sub-process (W*) + — > £ + £ + {W*)~ appearing in the diagram 
of Fig. 0. J ac and J^, are the quark (anti-quark) currents that couple in the t-channel to 
the virtual gauge bosons of the standard model (Fig. [I]) while J a ^ and J* d , are the incoming 
and outgoing currents of the qq' pair that couples in the s-channel to the W-bosons (Fig. [|): 



J^c = u(pc) r Pl u{p a ) J^ d = v(p b ) r Pl v{p d ) 
Ja, b = v(pb) r Pl u( Pa ) ( Jy = u( Pc ) r Pl v(p d ) 

with P L = (1- 75 )/2. 

With the following definitions of propagator factors: 

1/A = [{p a - Pc ) 2 - M 2 W ] [{p h - Pd f - M 2 W ] 
l/B = [{p a - p d f - M 2 W ] [(p b - p c f - M&] 



(14) 



l/A = [( Pa + p b ) 2 -M^ + iM w Y w ] [(p c + p d ) 2 -M^ + iM w Y w ] 



(15) 
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C = {p a + Pb ~ Pe) 2 ~M 2 N + iM N T N D = ( Pa +p b - p f ) 2 -M 2 N + iM N T 

the amplitudes (in the unitary gauge) are : 
(i) UiU, -► D k Di + £ + £ + 
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(19) 



where Ui denotes a positively charged quark (up-type) while Di denotes a negatively charged 
one (down- type). The quantities Vu i D- are the elements of the CKM mixing matrix. Of 
course the annihilation diagram of Fig. 2 comes in only in quark- ant i- quark scattering. 
In processes (i) and (Hi) the part of the amplitude depending on the factor B is due to 
the diagrams obtained upon exchanging the final state quarks. In the framework of the 
Lagrangian c.f. Eq. (|[), as discussed in section II, we have: 
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Due to the chiral nature of the couplings involved, the calculation is particularly simple if 
performed in the helicity basis |3(J. Within this approach the amplitudes are evaluated in 
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terms of scalar spinor products instead of scalar products of particle's momenta as in the 
usual method of squaring the amplitudes. In the massless approximation (all p 2 — for 
external particles) the only non zero scalar products between spinors of given helicity are 
the quantities: 



s(m,n) = s(p m ,p n ) = u + (p m )u^(p n ) , 
t(m,n) = t(p m ,p n ) = u-(p m )u + (p n ) , 



(21) 



which are given by: 



s(m,n) = -2^jE m E n G r , 



t(m,n) = +2\/E m E n R 



m J -"n M mni 



(22) 



with: 



G mn = cos(# m /2) sin(fl n /2) e +*(^-*»)/ 2 - sin(# m /2) cos(^/2) e -(<^-<M/ 2 , 
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(23) 



Consider for a moment the lepton line and its corresponding tensor T^ v projected out on 
the helicity basis: 



T^ x '\p e , Pf )=u x (p e ) 



l^Tiv . ivin 
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C D 

It is easily shown (use v\ (p) = —2X'j5U-x(p)) that only one of the four helicity combinations 
is non vanishing: 



T^' + \p e ,p f ) = T^'-\p e , Pf ) = 

T^' + \p e ,Pf) = T$'-\p e ,p f ) = u + (p e ) [2^ + *&] u _ {p f) 

In particular it follows that if p e = pf = p then C = D and: 

T l£~\p,P) = 2 -^r u + (p)u-(p) = s(p,p) = 0. 



(25) 



(26) 



Due to the structure of the coupling the two massless (identical) leptons have the same 
helicity and therefore by the Pauli exclusion principle if they have the same momentum 
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the amplitude must vanish. All amplitudes given below satisfy this property. The above 
remarks apply as well to the quark lines and in the massless approximation only one helicity 
amplitude is non zero. The following result is found: 
(i) UjUj -> D k Di + i + £ + 

~s(e,a)s(b,f) s(f,a)s(b,e) 



[n 
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(29) 



The above simple analytic form of the amplitudes is also quite easy to implement in a code 
for numerical applications, since the quantities s(pi,Pj) and t(pi,pj) are just functions of the 
energies and angles of the particle's momenta, c.f. Equations ([22] ) and (|23|) . 

In addition it is also possible to derive from the above expression of the amplitudes the 
average squared matrix elements for the various processes in terms of scalar products of 
particle's four momenta. One finds using \s(pi,pj)\ 2 = \t(pi,Pj)\ 2 = IpvPf 

(i)u l u J ^D k D l + e + e + 



Y,\M\ 2 = 128 K 2 \V* iDk V* ]Di \ 2 p c - Pd \p a -p e p b -p f 



(ii)U i D j ^D k U l +£+£+ 
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where the following definitions have been used: 



*i=(£ + Si) (s + £#)* & = ^a^W(%^ 



x 2 
x 3 



c + £25] (^d + £2^ J 6 = v UxD v UlDk /{y UiDk v UlD] ) 

C ~TS3 F \ D + ^3 1 



£3 - v UlDi v UkD j(v UkDi v UlDj 



(33) 



L(Pl,P2,P3,P4) = PVP2P3'P<1- PVP3PVP4 + PVP4P2-P3, 

£(Pl,P2,P3,Pi) = £^ XP (Pi)u(P2)v(P3)\(P4)p 



(34) 



with e 0123 = +1. The above formulas (|30|) , (|3l"r) , (32) have been written down in all generality 
allowing for possibility that the matrix elements of the CKM matrix might be complex. 
Except for s-channel contributions, such squared matrix elements were already given in 
fill in the framework of the same model as that here used. However, the formulas found 
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16| disagree from ours in various aspect, and in the fact that they seem to have been 
obtained as if a (1 + 7 5 )/2 were used for the coupling between leptons, Majorana neutrinos 
and W~ , instead of the usual (1 — 7 5 )/2 for left-handed coupling (formula (2) in [|T^] should 
then contain a (1 + 7 5 )/2 factor). The numerical applications reported here concern only 
quarks of the first and second family and for those the elements of the CKM matrix can be 
approximated (Wolfenstein parameterization ||31|| ) to a high degree of accuracy to be real. 
Accordingly ^1,^2,^3 will be real and QXi = QX 3 = 0. Only X2 will retain an imaginary 
part due to the complex resonant propagator factors. 
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The decay width of the heavy neutrino, T^, which appears in the factors C and D 
is of course a quantity which depends on the parameters of the particular model that is 
being considered here, B e ^ and Mat. It receives contributions from the two decay channels 
Ni — > W + P. (charged current) N% — » Z + vg (neutral current). In the present model there 
are no other decay channels open. Some authors however [|17| , |18H include also possible Higgs 
interactions of the heavy neutrino states. Doing so will open up additional decay modes 
depending on the mass relation between M^ and Mh- As both of these parameters are still 
unknown the present authors prefer not to consider this possibility. Hence the values used 
for Tjv in the present estimates of the resonant cross-sections should be regarded as lower 
bounds for the heavy neutrino total width. As the total decay width appears in the resonant 
propagator factors, it is a quite important parameter as regards numerical applications. 
Thus in appendix A are reported the formulas used to evaluate the partial and total decay 
widths, while in Table I the corresponding numerical values of the total decay width Tn 
are reported for different values of the neutrino mass M^. It can be inferred that the total 
decay width changes quite drastically going from T^/Mn = 0.0011% at M^ = 100 GeV to 
Tn/Mm = 0.43 % at M^ = 1000 GeV. In addition it is to be remarked that the authors 
of |7j quote even smaller values of the total decay width (see Table I) which appear to 
be unsubstantiated within the present model. Using the interaction of Eq. (§) one readily 
evaluates in the unitary gauge the partial widths (assume M N > 100 GeV) T(Ni — > W + I) 
and T(Ni — > Z + v{) and hence the total decay width. In the appendix it is shown explicitly 
that the values of the width quoted by the authors of ref. |7| can be reproduced by calculating 
in the 't Hooft gauge but not including the ghosts (Goldstone bosons) diagram. 

In reference || in order to keep the numerical computations of cross-sections reasonably 
simple, a constant value of T^ was assumed. Here this is not possible as the width of the 
heavy neutrino varies so strongly between low (« 100 GeV) and high (« 1000 GeV) masses, 
and one is obliged to keep the functional dependence of the total width on the neutrino 
mass Mjy. This makes the contribution of the resonant diagrams strongly dependent on the 
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heavy neutrino mass M^. 



IV. THE THIRD ORDER PROCESS pp 



-t + + W 



This process has been the object of where the square of the amplitude and the 
subsequent phase-space integration were calculated by means of the CompHEP package ||. 
The amplitude of the third order process is: 

1-75 



M = JC'A'v( Pb )Y 



u(p a )T^[eUQ)r 



v:,b*m a 



£' = ( ~7= j V ud a eN 1V1 N 

1/2' = [{p a + p b f -M 2 V + iM w T w ] 
C' = (p e + Q) 2 -M 2 N + iM N T N , 



D' = (p f + Q? ~M 2 N + iM N Y 



N- 



(35) 



with T^y defined in Eq. (p0|). Clearly relative to the annihilation amplitude of the process 
u d — y £ + £ + + ud, given in Eq. fll8l) , the outgoing quark's current [JY cd \\* is replaced by 
the on-shell W gauge boson polarization vector [cq\(<3)]* and the constants A' and K! are 
accordingly redefined. 

Here a formula is given which was derived by means of the standard method of squaring 
the amplitude and summing over polarizations of initial and final particles thus reducing 
the calculation to the product of traces of strings of Dirac matrices. 
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(36) 



where Q is the W boson four-momentum. It should also be remarked that the above result 
has been verified against that of the annihilation part of Eq. fl28|) by taking the small width 
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approximation on the latter, factorizing the phase-space and integrating out that of the two 
quarks from the on-shell W. 

However the formula given by the authors of ref. |7|] (c.f in the appendix) appears not to 
take into account the fact that the amplitude is the sum of two diagrams antisymmetrized 
with respect to the variables of the two identical leptons in the final state. The two diagrams 
have different heavy Majorana neutrino propagator factors (c.f C and D' in Eq. fl35f) ). 
Moreover all factors explicitly dependent on (M^) -1 and coming from the sum over the 
real W polarizations Ea^W^C 1 ) = ~Vv + {k tJL k u )/M^ v \ are missing. They claim that 
their result is the output of the CompHEP software package || which allowed them to 
perform symbolic computation of matrix elements. In order to understand the reason of 
this discrepancy the present authors decided to verify again the result given in Eq. ^ by 
using the same CompHEP software package used by the authors of ref. f?|. It appears that 
in CompHEP non standard models have to be defined, in the sense that new particles and 
new interaction vertices have to be defined by the user. CompHEP can perform calculations 
both in the unitary gauge and in the 't Hooft gauge. Defining the heavy neutrino and its 
interactions as in Eq. (|3|) and working in the Unitary gauge it has been verified that the 
CompHEP output reproduces that reported in Eq. ( |5B| ) up to the factor containing the 
totally antisymmetric Levi-Civita tensor: this is because the program directly replaces the 
squared propagator for the resonant particles. However, adding this term to the CompHEP 
fortran output of the matrix elements, it is found numerically irrelevant as it is proportional 
to the width of the heavy neutrino N. It is thus established that Eq. (|36|) agrees with 
the CompHEP output in the unitary gauge up to terms proportional to the width of the 
resonant particle (neutrino). 

Working with CompHEP in the 't Hooft gauge requires some care as in this gauge 
also the couplings of the W ± and Z Goldstone bosons (ghosts) have to be defined by the 
user. These are necessary to reproduce the correct sum over polarizations for W^ and Z. 
With the inclusion of the ghost vertices the calculation in the 't Hooft gauge reproduces the 
result of the Unitary gauge (as it should) which agrees with that of Eq. fl36D . Doing so the 
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present authors realized that the result of ref. |7j [c.f. Eq. (Al)] is indeed obtained with the 
CompHEP package working in the 't Hooft gauge if one does not defines the ghost vertices. 
The discrepancy of the result of this report, Eq. ([36|), with that of Eq.(Al) of ref. |7j is thus 
resolved: the formula given in ref. is just one of the terms when working in the 't Hooft 
gauge; the ghost contribution needs to be added to get the complete result which agrees 
with that reported here in Eq. (|36D which is obtained in the unitary gauge. 

In the appendix formulas are presented which illustrate this point explicitly for the sim- 
pler case of the heavy neutrino width. The calculation without including the ghost diagrams 
amounts to use [J2\ e >>(,k) e \(.k) = ~V^} f° r the sum over the gauge boson polarizations and 
the width thus obtained reproduces the anomalously small values reported in ref. for 
M N = 100, (800) GeV. 

The total cross-section, summing up the contributions of the quarks of the first and 
second generations, is given by: 



a = — dx a dx b ^2 [fu(x a )f£,(x b ) + x a ^ x b )] 

S J UD 



x I 4 £ \ M ^ UD "> e+e+ ^ 4LIPS(v/ ^ e ^ g) , (37) 
spins 

where the numerical factors account for average on initial spins and colors. All numerical 
results presented in the next section will be at the level of a parton Monte Carlo, while 
aspects like jet hadronization and detector simulations are beyond the scope of this work. 

V. DISCUSSION AND RESULTS 

The previous observation clearly does not affect the W-W fusion matrix elements because 
in the t-channel there is no finite width effect and the gauge boson propagator's numerator 
is simply t]^ u also in unitary gauge as the limit of massless external particles is assumed. As 
stated, correctly, in ref. these diagrams give rise to a very small cross-section, a ~ 10~ 5 — 
10~ 4 fb in the range of neutrino masses of interest here 100 GeV < Mn < 1000 GeV. At 
the LHC assuming an integrated luminosity over one year of L = 100 fb -1 one could define 
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a minimal cross-section cq = 10 -2 fb that would give 1 event per year. Clearly, smaller 
cross-section would be impossible to measure (with statistical significance) on a one year 
run. Thus, for the model studied in this work, the WW-fusion mechanism is negligible and 
not measurable. This was checked : (i) with authors' defined fortran code [|| based on 



Eqs. ( PD|pi|J5^ ) and feeding them into the phase space integration routines provided within 
the CompHEP package; (ii) and directly with CompHEP performing a symbolic matrix 
element calculation followed by numerical (VEGAS) phase space integration. 

Therefore the following discussion will concentrate on the resonant production mecha- 
nism of Fig. H], with a particular emphasis on the third order process and the comparison with 
the result reported in |7[]. The following numerical results were obtained with CompHEP, 
using the same kinematic cuts of ref. [0]: p^ pton > 5 GeV, p^ > 15 GeV and \r]\ < 2.5 
for the final state particles. The parton distributions functions (PDF) used are the set of 
MRS[G] p2| . As shown explicitly in these type of calculations are not very sensitive to 
the particular choice of PDF as long as one restrict itself to the more recent sets. 

Notwithstanding the differences highlighted in the previous section regarding the analytic 
formula of the matrix element and of the total width, the numerical value of the total cross- 
section is within rs 20% of that reported in |7| as can be seen in Fig. ^: this can be easily 
understood as the cross-section could be written, to a good approximation, as a product of 
the production cross-section a(pp-^£ + N) (independent from T^), times the branching ratio 
T(N— >£ + W~)/Tn that is numerically the same although the formulas are different. 

Differences are however found on angular distributions of decay products (leptons and/or 
jets), that may be very important for the detection of the signal. Notice that from theoretical 
arguments |33| it can be stated that a forward-backward symmetry should show up in 
processes of production of heavy Majorana through CP conserving interactions like that 
here considered. Figs. |]&[| show normalized angular distributions of the outgoing lepton(s) 
and gauge boson. Let Qi (Qw) be the angle formed by pg (pw) with the positive direction 
of the collision axis, in the center of mass system (c.m.s.) of the hadronic (pp) reaction 
(lab frame). Figs. f|&|5] show the expected forward-backward symmetric behavior and in 



particular that the final state particles are emitted preferentially parallel or anti-parallel to 
the collision axis. Fig. ^| shows the normalized distribution in the opening angle of the two 
like sign leptons. Here a difference is found relative to the numerical results reported in 
ref . |7j . In Fig. |6| one can observe that for relatively low neutrino masses the two leptons are 
emitted preferentially parallel to each other while at high neutrino masses the distribution 
tends to become flatter. 

The opposite behaviour is instead predicted in ref. |7| (see their Fig. 5). This can again 
be explained assuming the interpretaion given above of the calculation reported in [0]. The 
absence of the ghost particles amounts to the absence of the longitudinal polarization of the 
W-gauge boson (J^ = 0): so the W gauge boson would have only the two helicity states 
A = ±1. Suppose a head-on collision between a u and d with p u along the positive direction 
of the z axis. As u has A = — \ and d has A = +|, the initial total angular momentum is 
J z = — 1. To conserve angular momentum along the z axis, the only allowed configuration in 
the final state, composed of the two identical leptons and an on-shell W, is with a W-gauge 
boson with A = — 1 and the two positrons emitted with opposite momentum (antiparallel), 
since the two leptons have the same helicity, to give J*f = 0. If we start with d u we have 
an initial J z = +1, so the final W has A = +1 and the two leptons are again emitted 
antiparallel. However in reality the W-gauge boson has also the longitudinal polarization 
and the JY = possibility is also present: when this is so, the two leptons must be parallel 
to give J l z l = ±1. This additional helicity amplitude prevents the angular distribution to 
drop as 9 U — > 0. 

The other distributions discussed in ref. (c.f. transverse W-mass, transverse W- 
momentum and invariant mass Mew, Muw) are confirmed by the present calculations and 
so are not repeated in this report. Fig. |7] shows the angular distribution with respect to the 
opening angle of the jets in the fourth order process. It is worth noticing that the W-angular- 
distribution, Fig. [5], is more strongly peaked along the collision axis than that of the leptons, 
Fig. [|. Thus the jets that originate from its decay are expected to be emitted preferentially 
along the collision axis and parallel to each other. This is indeed confirmed by the parton 
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level simulation as shown in Fig. [I]. If one takes into account that in the hadronization 
process more than two jets may be present, it is clear that the mass reconstruction of the 
W may not be so easy from the experimental point of view. 

The peculiar characteristics of the signal are: two like sign leptons with the same angular 
and transverse momentum distributions, no missing energy, and at least two jets with an 
invariant mass distribution peaked at the W mass while the invariant mass distributions in 
the jj£ system will be peaked at the heavy neutrino mass, as shown in @]. Clearly no SM 
process have these features but some have final states with accompanying neutrinos (missing 
energy) which are a source of potential background to the LSD + 2jets signal. All previ- 
ously cited literature considered only heavy quark (it, bb) production and subsequent decay 
as a major source. However recent simulations |34[ have recognized that more dangerous 



background at LHC can be single top production and gauge boson pair production (WZ 
and ZZ) followed by leptonic decays. After all cuts (for details see ref. p4|| ) the authors find 
a background of 4.9 ±1.6 events with an integrated luminosity of 10 fb _1 that implies a 
total cross a ~ 0.5 fb. This would restrict the observability of the signal studied here in the 
neutrino mass range 100 — 250 GeV. However the kinematic cuts used here are softer than 
those used for the SM background suppression in ref. |34| , so it is interesting to understand if 



is 



the conclusions above can be significantly affected by changing the kinematic cuts. In |34 
indicated that at least a cut of 40 GeV on p^P on is necessary to suppress the SM background 
(also an isolation cut is needed but this will not affect the signal significantly). In Fig. |3| 
the dashed curve shows the effect on the signal of increasing the cut on the p^ p on up to 40 
GeV: the total cross-section is lowered by a factor of up to an order of magnitude in the 
region of low neutrino masses, but it is practically unchanged in the high mass region. The 
drastic change in shape found at small masses can be qualitatively understood by noticing 
that when M^ is small a larger part of the energy goes into My/ mass leaving not much 
kinetic energy for the outgoing particles, so a high cut on pt naturally lowers the cross- 
section. This is in agreement with the lepton transverse momentum distributions (which 
are not reported here as they are the same of those reported in ref. 0): these are peaked 
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at Pt P ° n ~ Mjv/2. Clearly a cut of 40 GeV on p^P on cuts away most of the events when 
M N ~ 100 GeV, thus causing a drastic change in a tot . However it should be remarked that 
also within these strong cuts there still remain a measurable number of events. 

VI. CONCLUSIONS 

This paper has addressed the lepton number violating processes (i) pp-*£ + £ + + 2jets 
and (ii) pp — ► £ + £ + W~ within the scenario of models where heavy Majorana neutrinos mix 
with light leptons. In particular regarding the fourth order process (i) explicit analytic 
expressions of the helicity amplitudes are provided for the two competing mechanisms: the 
WW-fusion and the resonant neutrino production. These helicity amplitudes might also be 
useful for more realistic (detector dedicated) Monte Carlo simulations and/or within the 
scenario of other models beyond the standard electroweak theory. With respect to process 
(ii) a detailed study, both analytical and numerical, with a thorough comparison to the 
work of ref. JF|, has been performed finding some differences in the analytical formula of 
the matrix element which have been explained by noticing that the formula given in ref. |7j 
does not include the ghost contribution needed in the Feynman-'t Hooft gauge. However, as 
important as this might be, the main conclusion of ref. |7|] (the size of a to t) is only affected by 
an error of ~ 20%, though relevant differences are found on some of the angular distributions. 

In particular the distribution in the opening angle (0«) of the two leptons (Fig. ||) is 
found to be significantly different from that reported in |7j . 

According to the study presented here in the framework of the mixing model of heavy 
left-handed neutrinos there is a mass range of up to M^ ~ 800 GeV where the cross-section 
would be in principle observable ( 1 event per year assuming an integrated luminosity of 100 
fb _1 ). This window in the heavy neutrino masses is drastically reduced down to Mn ~ 250 
GeV taking into account the recent SM background estimates quoted in ref. [|4]] (~ 50 ± 16 
events/year with an integrated luminosity of 100 fb _1 ), which appears to be dominated by 
the production of double gauge bosons and single top. 
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The observation of the LSD signal, in this model, is thus ultimately related to the size 
of the mixing angle and the possibility of suppressing the SM background. If a like sign di- 
lepton signal will be seen at LHC, then the question of determining the underlying physics 
scenario, that triggers the non SM signal, would immediately arise; more realistic, detector 
dedicated, Monte Carlo simulations would be needed and, of course, this eventuality would 
require additional theoretical work. It is the authors' belief that the present work would in 
any case be of importance in this direction. 
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APPENDIX: TOTAL DECAY WIDTH OF THE HEAVY NEUTRINO 

The expression for the partial decay widths of the heavy neutrino states as deduced from 

the interaction Lagrangian given in Eq. (|8|) are given by: 

n A/f 3 
T (N t -> Wt] = — — -f \B eN \ 2 $ (x w ) , (Al) 

16 s w Mfo 

n A/f 3 
T(N g - Zu t ) = — — -f \B eN \ 2 \B vN \ 2 $ (x z ) , (A2) 

16 s w Mfo 

with $(a;) = (1 — x) 2 (l+2a;), Xw,z — ^wz/^n anc l s w — sin 2 ®w-, where 9w is the Weinberg 
angle, a the fine structure constant and \B V ^\ 2 is determined from Eqs. (§). 
The total decay width is given by 



N 



2 T(N e -> Wt) + T(N e -+ Zu t ) (A3) 

The factor of two in Eq. ( |A3| ) arises because, due to its Majorana nature, the heavy neu- 
trino can decay into W~£ + and W + £~ with equal probability. In Table I, the corresponding 
numerical values are reported for the total decay width T^. Including Higgs interactions of 
the heavy Majorana states will increase the neutrino total decay width T^. 
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The formulas that reproduce the widths given in 0, as explained in the text, are obtained 
using 't Hooft-Feynman gauge but without including the Goldstone boson contribution, 
which amounts to use for the sum over the real W polarizations: J^ A e\(k)e* x (k) = —rj^, i.e. 
that corresponding to a massless gauge boson. One finds: 

T{N, - W£) = - M N \B eN \ 2 U(x w ) , (A4) 

o Sw 
01 

T(N e -> Zv t ) = M N \B eN \ 2 \B uN \ 2 n (x z ) , (A5) 

with n(x) = (1 — x) 2 and c\y = cos 2 9\y, which indeed reproduce the values of T N given by 
and reported in table I. 
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TABLES 
TABLE I. Total decay width of the (electron) heavy neutrino T N = T(N -»• Wl)+T(N -► Zv), 



obtained through the use of Eq. ( |AlD and Eq. ( jA2|) with 5 eA r = 0.0052. 



Mjv(GeV) 100 300 500 600 800 1000 

LAr(GeV) 0.11 x 10~ 2 0.138 x 10+° 0.65 x 10+° 1.13 x 10 +0 0.267 x 10 +1 0.52 x 10 +1 
LAr(GeV) a 0.650 x 10" 3 - - - 0.580 x 10" 1 



1 For comparison the values used by the authors of M are also reported. 
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FIGURES 

FIG. 1. W fusion mechanism for production of Like-Sign-Dileptons (LSD) in high energy 
hadronic collisions. These diagrams, including the appropriate exchanges, describe also quark 
quark scattering and anti-quark anti-quark scattering. 

FIG. 2. Production of LSD through quark anti-quark scattering. In addition to the diagram of 
Fig. 1 (virtual W fusion), one must consider the £ + Np production via quark anti-quark annihilation 
with subsequent hadronic decay of the heavy neutrino Ni — > £ + qq {£ = e, /U, r). In principle for qq 
this process interferes with the one in Fig. 1. 

FIG. 3. Signal cross-section for pp^£ + £ + W~ (£ = e): the solid line includes contributions 
from first and second generation quarks; the dot-dashed line includes only the ud contribution 
which is clearly the dominant (90%); the dashed line indicates the total cross-section when a 
higher cut on the lepton transverse momentum px > 40 GeV is applied. 

FIG. 4. Lepton polar angle (@i) distribution in the third order process pp — > e + e + W~ at LHC 
energy v5 = 14 TeV, for different values of the heavy neutrino mass M^. @e is the laboratory 
polar angle i.e. the angle formed by the lepton with the positive direction of the hadronic collision 
axis. 

FIG. 5. Distribution in the gauge boson W polar angle (@w) m the third order process 
pp — > e + e + W~ at LHC energy \/S = 14 TeV, for different values of the heavy neutrino mass Mjy. 
The same legend as in Fig || applies. The polar angle @\y is defined in the labratory system i.e. 
with respect to the hadronic collision axis. 
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FIG. 6. Distribution in the opening angle between the two identical leptons (©a), defined in 
the laboratory frame, in the third order process pp — ► e + e + W~ at LHC energy yS = 14 TeV. The 
legend of Figure || applies. 

FIG. 7. Fourth order process pp — ► e + e + 2jets. Distributions for the opening angle (in the 
laboratory frame) between the two quarks (jets) (®jj) coming from the W~ gauge boson decay. 
The different curves refer to different values of the heavy neutrino mass according to the legend of 
Figure [|. 
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